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1.0 SUMMARY 


Tills User’s Manual contains a complete description of the computer codes known 
. ,c AX.SVMMKTRIC D.m.SRR W.CT cod. or ADO rode. U lorlodooTuot o^;",v;o;:oo 
Which describe the formulation of the ADD code and comparisons of calculation with 

mT' Z :y ’“"■r.”"’ of fW cod.. 1 . doocllkod ir 

ttu first volume. The second volume contains a detailed description of the code 

Including the global structure of the code, list of FORTRAN variables, and descrip- 
ulU.T code which generates coordinate systems for arbitrary axi8>*mmetrlc ducts. 




2.0 INTRODUCTION 


Thin U«or'» Manual doHcrlbos tho computer cocIoh knovm collectively an tlu* 

AXl SYMMETRIC DIFFUSKR DUCT code or ADD code. This code was originally developed for 
NAJ>A Lewis Research Center under contract NAS3-15402. Important revisions, including’ 
the conformal mapping coordinate generator, were developed for the U.S. Army Air 
Mobility Research and Devt ^/ment Laboratory under Contract DAA%T02-73-C-0037 . Fur- 
ther developments and improvements to the ADD code were funded by United Technologies 
Research Center and Pratt & Wlutney Commercial Products Division. Additional improve- 
ments, including incorporation of a two equation model of turbulence and a compres- 
sible axisymmetric streamline curvature correction was funded under NASA contract 
NAS3-218S3. Finally a new coordinate generator which permits 180 deg turns in ‘ duct 
was developed under NASA contract DEN3-235. 

All the important features of the basic analysis contained in the ADD code have 
been published in the open literature. Tl\e accuracy and reliability of the code has 
been demonstrated by published comparisons of calculated results with experimental 
data. The basic analysis used in the ADD code was developed by Anderson (Ref. 1 and 

A description of the blade force calculation is given by Barber et.al.(Ref. 3). 
Tlu' code has been successfully applied to predicting the performance of the subsonic 
portion of mixed compression inlets by Bowditch (Ref. 4) and to predicting the pres- 
sure recovery of high Mach number diffusers by Povinelli (Ref* 5). Additional appli- 
cat ions of the ADD code have been to analyze swirling flow in a precombust ion diffuser 
and also the flow in an inlet with inlet guide vanes (see Barber et al. Ret. IV 
Finally flows in small gas turbine ducts have been analyzed by McLallin and Kofskey 
(Ref. 6). Tlie three turbulence models incorporated into the ADD code and also the 
compressible axisymmetric streamline curvature cc^rrections have bec'n evaluated by 
Anderson and Edwards (Ref. 7). Modified versions c^f the ADD code have been developed 
to analyze external fU^ws such as underexpanded hot supeirsc^nic jets expanding into 
cold subsonic mainstreams by Vatsa ot.al (Ref. 8) and alsc^ the high speed flc^w inter- 
act ivMi between a propeller and nacelle by Kgolf et al.(Ref. ^). 

This User's Manual has Deen organized into several sections for the convenience 
of the users. Section 3.0 contains a descript ic^n of the different versions of the 
ADD ccxle and a description of their special features. This sectic^n should assist 
the general user in determining if the ADD code is applicable to the problem at hand. 
Section 4.0 contains a detailed description of the c^peration of the ADD code, includ- 
ing a typical run stream for UNIVAC computers, input/output formats and sample prob- 
Itnns. In addition to operation of the code, this section contains a list of DIAG- 
NOSTICS, which are internal checks within the code to measure the progress of the 
calculation* If the code fails and prints a DIAGNOSTIC, this section should assist 
the user in determining the cause of the failure and suggest a nmedy. Sections 5.0 
through 7.0 are written for the special user who wishes to modify or upgrade the code 
tor a particular problem. This portion of the manual contains sections on the general 
structure of the code, definitions of COMMON block variables, and detailed descrip- 
tions ot each ol tiu' subroutines. Sections 10.0 through 12.0 contain a descrip- 
tion of the CODUCT code which is an alternate mesh generation code developed for NASA 
under Contract Dr.N3-23S (Ret. 10). 
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3.0 GENERAI- DESCRIPTION 


3.1 Versions of the ADD Code 

Four versions of the ADD code are currently in general use. The basic code is 
called the ADD code and is described in Ref. 1, 2 and 3. The version of the code 
called the PREMIX code was developed for NASA Lewis Research Center under contract 
NAS3“21269. This version of the code was developed to analyze the performance of 
premixing prevaporizing fuel air mixing passages. It consists of three codes; ADD, 
PTRAK, and VAPDIF. The PREMIX/ADD code differs from the basic ADD code only in in*” 
put/output which is in International Standards units. PREMIX/PTRAK code solves the 
problem of tracking vaporizing fuel droplets in a three dimensional flow field. 
PREMIX/VAPDIF solves the problem of discussion of fuel vapor into a moving air stream. 
These codes are described by Anderson et al. (Ref. 10) and applications of these 
codes to specific premixing passages are given by Anderson et al. (Ref. 11). The ADD/ 
JET code is a version of the ADD code which is modified to treat the expansion of a 
hot underexpanded supersonic jet into a cold subsonic free stream. This version of 
the ADD code differs from the basic ADD code in the use of a computational grid and 
model more suited to jet flows. In addition it uses streamline curvature 
obtained from a separate calculation. A description of the principal features of 
this code and a comparison of calculated results with experimental data is given by 
Vatsa et al. (Ref. 8). The PROPFAN version of the ADD code was developed for NASA 
Lewis Research Center under contract NAS3-20961. A description of this code and 
preliminary results are presented by Egolf et al. (Ref. 9). This version of the ADD 
code was developed to treat the high speed flow interaction between a propeller and 
nacelle and differs from the basic ADD code in the use of a propeller lifting line 
analysis which is used in place of a compressor cascade analysis. 
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3.2 General Features of ADD Code 


Program Languag e 

The Axlsynunetric Diffuser Duct (ADD) code source program is written in FORTRAN V 
computer language for use on a UNIVAC 1100/81A computer. Some machine specific 
language, such as PARAMETER and INCLUDE FORTRAN statements is used. However, these 
statements may be replaced easily by eqivalent code for use on other machines. 
Successful conversion o^ the code to both IBM and CDC computers has been made 
and these versions of the code are available. The ADD code makes use of a UNIVAC 
routine NTRAN which stores and retrieves large data blocks on disc files; however, 
the ADD code Is organized so that NTRAN is easily replaced by the equivalent 
F0RTRAN DEFINE FILE. Finally, it should be noted that the ADD code makes use of 
least squares spline fitting and smoothing subroutines provided by IMSL, Inc. 
which are available at all major computer centers. 

Types of Fluids 

The ADD code can treat any compressible fluid with constant thermodynamic 
properties for the gas constant R and the specific heats Cp and Cv. The molecular 
viscosity, which is temperature dependent, is estimated using Sutherland's law; 
the molecular thermal conductivity is calculated using a constant value for 
Prandtl number. The viscosity of the fluid at standard conditions and Prandtl 
nvimber are input parameters. If these properties are not specified in the input 
data, the ADD code uses the properties of air at standard conditions. 

Types of Flow Treated 

The ADD code may be used to treat any subsonic compressible laminar or turbulent 
swirling flow in axisymmetric ducts or nonswirling flow in two-dimensional ducts. 

The duct shape may be annular or two-dimensional with both inner and outer walls; 
or, it may be an axisymmetric duct with only an outer wall. Subsonic flows have been 
calculated successfully up to choked conditions. The mixing of hot and cold flows 
in a duct have also been calculated successfully. The code, however, cannot 
calculate flows with significant regions of separated or reverse flow. 

Duct Geometry Options (I0PT3) 

The flow through any axisymmetric or two-dimensional duct can be calculated. 
Ducts with sharp discontinuities in flow area, which produce flow separation cannot 
be calculated. 

For convenience, provision is made in the code to analyze flows in straight 
annular ducts (I0PT3*1) or in straight wall, annular diffusers (10PT3*3) using only 
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a few input parametern. For ducts of arbitrary shape (I(HPT3»2), the coordinates 
(radii) of the inner and outer walls are specified at JLPTS equally-spaced axial 
stations. To cssure that the curve representinq the duct contour has continuous 
first and second derivatives, a least-squares spline fitting, smoothing and Inter- 
polation procedure is Included in the code. This procedure is used whenever the 
number of streamwisc stations (JL) is not equal to JLPTS. When the new C0l)lifT code 
is used to generate coordinates 10PT3"4 

The spec IflcMt ion of the duct geometry must include a straight, annular inlet 
section whose length is at least equal to its height. Two-dimensional ducts are 
treated as annular ducts In which the height of the duct Is small compared to the 
radius of the duct. Numerical experiments have shown that, if the height of the 
duct is less than 1/100 of the duct radius, the flow is essentially two-dimensional 
to an accuracy of three decimal places. 

Inlet Flow Options (10PT1) 

Any arbitrary inlet flow conditions may be specified which is consistent with 
the equations of motion and the turbulence model. Two types of input data are 
required; (1) specification of the Inviscid free stream and core flow conditions, 
and (2) specification of the laminar or turbulent boundary layer flow parameters. 

With (I0PT1*3, A, 9), the flow is assumed to be turbulent and with (I0PT1«7, 8) 
the flow is assumed to be laminar. With I0PT1»3 or 7, the core flow is calculated 
assuming that the stagnation pressure and stagnation temperature is constant 
across the duct. The input Mach number and swirl angle determine the velocities 
and weight flow, and the static pressure is determined by solving the conservation 
equation for radial momentum. When 10PT1*>4 or 8, the inlet core flow is determined 
by specifying KLL data points for fractional distance Y, stagnation pressure PfCY), 
static pressure P(Y), swirl angle a(Y), and stagnation temperature TjfY). For 
10PT1*9, the core flow is determined by specifying KLL data points for fractional 
distance Y, streamwise velocity Ug(Y), stagnation pressure PfCY), swirl velocity 
and stagnation temperature Tj(Y). Isentroplc flow relations and radial 
momentum conservation equations are used to determine the remaining variables. In 
addition, when 10PTL“A, 8 or 9, the corresponding exit flow data must be provided. 

If the exit plane data is not available, the inlet plane data may be repeated. 

The boundary layer velocity and temperature profiles are constructed from known 
analytic solutions using the boundary layer displacement thickness <S* and a power law 
(1/n) velocity profile. For laminar boundary layers (I0PT1-7, 8) a Balsius profile 
Is assumed. For turbulent flows (10PT1-3, A, 9), Cole's boundary layer profile Is 
used with the shape parameter determined for 1/n. 

In many flow situations, it is often more convenient to specify the weight flow 
rather than velocity or Mach number. For these situations, the user may specify 
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^ "kI' »• •">» Pr«.ure profile 1. auto,»atlc»ll>. 

odjuotod to obtain the raqulrad valght flou with the other Input varlahlee held 

llxed* 


n*^ **^°r^*^ *^^‘**^ Initial i-lane condltlona muat aatlsfy tho Xawa of 

motion and be compatible with the turbulence model. Therefore, the ADD code makea 
many checks on the input data to Insure satisfactory starting conditions. As an 
example, the Initial plane data Is checked to determine If the radial momentum 
conse^ation equation is satisfied. If it is not satisfied, the input static pressure 
profile is replaced by the static pressure calculated from the radial momentum 

initial" ® DIAGN0STIC message is printed. The weight flow calculated from the 
initial plane data is checked to see if it is greater than the choked-flow value. 

If it is greater, the calculation stops and the value of the choked weight flow 
« Checks are made to Insure that the boundary layer profile can be 

L thf 11 *'*“*?" necessary adjustments are made automatically 

and the calculation continues. In all cases where adjustments to the input data 

are made and the calculation continues, a D1AGN0STIC message is printed. When 
no adjustment is possible or when the flow situation is physically impossible, the 
calculation stops and the user is notified with a DIAGN0ST1C message. A list of 
these DIAGN0STIC messages is given in Section 4.4 

Grid Selection 


The user may determine the calculation grid using input parameters or the 
grid may be determined automatically. In either case the user must specify the 
number of streamwise stations (JL) and the number of streamlines (KL) . Experience 
has shown that a 50 x 50 mesh is suitable for most problems. Default options 
exist for both the distribution of mesh points In the cross flow direction as 
specified by the mesh distortion parameters DDS and the streamwise step size 
para^ter ^S. In selecting the mesh distortion parameters DDS, numerical accuracy 
requires that a sufficient number of mesh points exist in the turbulent sublayer. 

In practice, the first mesh point from the wall should be at Y+ » 1.0 and at least 
0 mesh points should be in the boundary layer. Since these criteria depend on 
both the flow Reynolds number and wall friction coefficient, they are not con- 
venient for the user to calculate a-priorl. Therefore, if DDS is not specified 
in the input data, a value for DDS is calculated using an algorithm which produces 
good results for most cases. The value for the streamwise step size parameters 
ros depends on the boundary layer thickness and race of growth of the boundary 
layer. If KDS is not specified, the code selects a value for KDS between each 

streamwise station using an algorithm which produces satisfactory results for 
most cases* 


Print Options (10PT4) 

The frequency and quantity of output are controlled by the print option 
I0PT4. If 10PT4 >0, the output consists of the mean flow variables including 
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fitroamwlfiii velocity Ug, CanKoncial velocity U., awlrl angle atagnatlon 
preanure Pj, atagnatlon temperature Tx» and Mach number M at each atrearowlae 
atatlon for JL atatlona; thla printout occurs every I0PTA*-*' atatlon. If 
I0PTA ^ -1, additional Information la printed Including the effective 
turbulent vlacoalty and thermal conductivity, the boundary layer aolutlon In 
universal coordinates U'^(Y^) , and the turbulent kinetic energy distribution; this 
Information la printed every I0PTA^** station. 

Diagnos tics 

The ADD code makes numerous checks during the progress of the calculation. 

If the program Is able to remedy a detected problem, a DIAGN0STIC Is printed and 
the calculation continues. If a fatal error Is detected, the calculation stops 
and a DIAGN0STIC notifies the user about the nature and location of the error. 

A complete list of D1AGN0ST1CS is given in Section A. A. 

Coo rdin ate O ption ( 10PT9) 

The calculation of the coordinate system may be stored on a data file and 
retrieved for use in subsequent cases. If I0PT9«1, both the coordinates and the 
viscous flowfield are calculated. If I0PT9*2, the coordinate system is calculated 
and stored on file NINE and the calculation stops. If I0PT5»3, the coordinates 
stored on file NINE are recalled and the viscous flowfield is calculated. This 
feature is particularly useful when the user wishes to calculate several flows 
using the same duct geometry. If CODUCT is used I0PT9*3. 

Data Files 


A list of data files and storage requirements are given on Table 1, Section 5.5. 
The ADD code or C0DUCT code generates two coordinate files. File NINETEEN is a 
coordinate file with a uniform mesh, and File NINE is a coordinate file with a mesh 
distorted to provide grid resolution in the boundary layers. In addition, the in- 
viscid flow field solution is stored on File IWENTYTW0 and the viscous solution is 
stored on file EIGHT. It is recommended that these files be registered and cata- 
logued so that the data may be stored permanently over a period of several weeks. 
Proper use of these files allows the user increased flexibility in solving problems. 

Start /Stop Options 

A flow calculation may be started at coordinate station J-I0PT15 and it may be 
terminated at coordinate system station J-I0PT16. If I0PT15 is not specified, it 
is assigned a value 10PT15*!; if I0PT16 is not specified, it is assigned a value 
I0PT16*JL. The calculation of the flowfield may be continued (or restarted) at the 
JM coordinate station by specifying 10PT17»JM only if in the preceding calculation 
10PT1A > 0. 
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T urbultmce Models (10PT12) 


The MID code is provided with four optlonnl turbulence models described In 
Ref. 7. For I0PT12'*O, I, 2 al^obrslc turbulence models nre used based on 
Prnndtl’s mixing length theory. For I0PT12«3, « two equation model of turbulence 
In used. Op tion I0PT12"O usen n, turbulence model which Is well estnbllshed for 
equilibrium turbulent flowftelds and Is therefore recommended for all calculations. 
The other options (I0PT12'"1, 2, 3) are oporatlonnl but these models liave b<*en 
applied to only a few flowflcld sltuotlona; the use of these models Is not 
recommended at the present time. 

Bl ade Force Options (1<I>PT2). (10PT5). (I0PT1(^ 

Struts, inlet guide vanes, stators, and rotors are modeled in the ADD code 
as a-priori body forces. Three options exist in the code for calculating these 
forces. If measurements of stagnation pressure P«y, swirl angle a, and stagnation 
temperature T^ are available, the blade forces can be calculated from blade element 
theory by settlag 1J>PT2=1. If I0PT5»2, the program uses the inlet/cxit flow data 
for I0PT1*4. It 10PT5«1, separate data must be loaded for the blade force 
calculation. If I0PT2>3, the blade force is calculated from the flow conditions 
and blade geometry using blade element theory and empirical cascade correlations. 

If I0PT2*4, the blade force is calculated using the distributions of exit air 
angle 02 (Y) and loss coefficient ZgfY). 

I0PT1O determines whether the blade is stationary (I0PT1O=O, stator) or 
rotating (I0PT1O»1, rotor). 

Splitie i tt ing Option JL f* JLPTS 

Man; ..ct contours can only be obtained by measuring coordinates from an 
engineering drawing. Since the ADD code requires curvature (i.e., second 
derivatives) these measured coordinates must be very accurate. In general practice 
this accuracy is not possible so therefore, a general spline fitting, smoothing, 
and interpolation routine is used. This subroutine makes use of a standard IMSL 
routine ICSVKU which is a spline fitting routine which optimizes the location of 
the knots or nodal points. The wall contour is numerically differentiated to 
obtain second derivatives. A spline Is fitted to the second derivative and 
integrated analytically. Thus the wall contour is continuous up to the fifth 
derivative. This option is used when the number of output data points JL does not 
equal the number of input data points JLPTS. 

Streamwlse Curvature Correction I0PT7 

When I0PT7“O, the ADD code uses the streamline curvature data stored on file 
NINE. When 10PT7*1, the ADD code calculates the compressible axlsymmetrlc potential 
flow solution, the corresponding streamline curvatures and stores the results on 
file NINE. Subsequent calculations can then be made with I0PT7=O. 
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A.O OPERATION OF ADD CODE 


A.l RunstrPrtm for ADD Cod*^ 

Tho following runHtroam In nufflclont to execute the ADD 
1100/BlA computer using the Exec. B operating nyNtern. 


0ASG,A 

RIGHT., D/O/TRK/ 500000 

OUSE 

8., EIGHT. 

(MSC.A 

NlNE.,D/O/TRK/500000 

«>USE 

9., NINE. 

C^ASG.T 

10,D/0/TRK/6000 

@ASG,T 

11,D/0/TRK/6000 

0ASG.T 

12,D/0/TRK/6000 

0ASG.A 

NINETEEN. ,D/0/TRK/500000 

0USE 

19., NINETEEN 

0ASG,A 

TWENTYTW0. .D/O/TRK/500000 

0USE 

22. .TWENTY TW0 

@ASG»T 

23.D/O/TRK/15000 

0ASG,T 

24.D/0/TRK/15000 

0ASG,T 

25,D/0/TRK/50000 

0XQT 

MAPADD 

Input Data 


0FIN 



ORIGINAL PAGE |s 
OF POOR QUALITY 


code on the UN 1 VAC 
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4.2 Input Format for ADD Code 

The input format for the ADD code la described on the input data coding forms 
which follow. These coding forms are organized with one form per input data card. 
Each form contains the names of the variables, the format, and a description of the 
data. The input option card controls the data that must be read. Since not all 
cards are read, the user should make certain that the input data agrees with the 
input options. 

In general the input data is read as follows: 

Card 1 Title Card 

Card 2 Option Card 

Card 3 Mesh Parameter Card 

*Card 4 Duct Geometry Card 

+ data as required by I0PT3 
*Card 5 Inlet Flow Card 

+ data as required by I0PT1 
Card 6 Force Data Card (If I0PT2 # 0) 

+ data as required by I0PT2, I0PT5, I0PT1O 
*Card 7 Reference Card 

Card 8 Slot Flow Data Card 

+ data 

Card 9 Wall bleed data card 
*Card 10 Interpolated output data card 


* NOTE: Blank cards must be loaded when options are not used. 

See detailed writeup. 
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I0PT1 


(FL0WIN Option) 


I0PT2 


I0PT3 


10PTA 


■ 3 Inlet flow is computed by specifying data on Card 5 (turbulent flow). 

- 4 Inlet and exit flow profiles are read from 2*KLL data cards 

following Card 5. Input, fractional distance Y, stagnation 
pressure Px<Y), static pressure P(Y), swirl angle a(Y), and 
stagnation temperature T^(Y) (turbulent flow). 

- 7 Inlet flow is computed by specifying data on Card 5 (laminar flow). 

• 8 Same as 4 but for laminar flow 

■ 9 Same as 4 but: Input fractional distance Y, static pressure 

Pl>(Y), streamwise velocity Ug(Y), swirl velocity U^(Y) and 
stagnation temperature T^(Y) (turbulent flow). 

(FORCE Option) 

• 0 No blade force 

• 1 Calculate blade force from upstream/downstream flow data; input 

fractional distance Y, stagnation pressure Pj(Y) , static pres- 
sure P(Y), swirl angle a(Y), and stagnation temperature T.j>(Y) 

• 2 Not available 

• 3 Calculate blade force from cascade correlations 

“ 4 Calculate blade force from fractional distance Y, exit flow 
swirl angle oofY), and loss coefficient Zg(Y) 

(GDUCT Option) Information follows Card 2 

• 1 Calculate a straight, annular duct 

“ 2 Arbitrary duct with evenly spaced axial stations 

- 3 Calculate a straight wall annular diffuser 
" 4 Coordinates stored on data file 

• 5 Arbitrary duct with arbitrary axial stations 
(PRINT Option) 

Print solution every I0PT4 station. For example, if 10PT4 ■ 3, every 
third station will be printed. Tf I0PT4 - -1, the code provides 
an extended printout; this extended printout includes information 
about tile boundary layer profiles and the turbulence model. 


i-n 


I0PT5 


(STRUT INPUT Option) 


I0PT6 


I0PT7 


I0PT8 


I0PT9 


Strut input data (If I0PT2 - 1) used to calculate strut forces 
from experimental data measured upstream and downstream of strut* 

■ 1 The upstream and downstream strut data cards are identical 
to the Inlet and exit flow cards and are not read. 

• 2 Read in required profiles. 

(STRUT Thickness Effects) 

0 Include strut forces plus thickness effects 

1 Include strut thickness effects only. 

(Axisymmetric Compressible Streamline Curvature Corections) 

0 “ No curvature correction 

1 - Curvature correction 
(WBLEED Option) 

» 0 No Bleed 
« 1 Bleed OD wall 
<■ 2 Bleed ID wall 

> 3 Bleed OD and ID wall 
(COORDINATE Option) 

> 0 Make an approximate calculation for both streamlines and 

potential lines— do not save flowfleld on disk. Used only 
for I0PT3-1 

■ 1 Make exact calculation of streamlines and potential lines — 

store results on logical unit 9 and complete viscous flow 
calculation. 

• 2 Same as 1 but terminate calculation after coordinate calculations 
are completed. 

■ 3 Read geometry from logical unit 9 and use in viscous flow 

calculation. 
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I0PT1O 


(R0T0R Option) 


I0PT11 


I0PT12 


I0PT13 


* 10PT14 


10PT15 


10PTlb 


I0PT17 


0 ■ Stator 

1 ■ Rotor 
(FL0W Option) 

• 0 Interal flow. 

■ 1 External flow. 

(TURBULENCE Option) 

• 0 Use two-layer turbulence model. 

* 1 Use two-layer turbulence model with low Reynolds number 

correction, (not tested) 

“ 2 Use two-layer turbulence model with streamline curvature correction. 

* 3 Use two equation turbulence model (applicable to flows in annular 

diffusers only; i.e., diffusers with both inner and outer walls). 

(SL0T Option) 

■ 0 No slot cooling. 

« 1 Slot cooling. 

(GL0BAL Option) 

•• 0 Global iterations not used. 

> 1 Global Iterations used - backward differencing for streamwise 
velocity derivatives in vicinity of separation. (See Ref. 7) 

(JFIRST Option) 

Start flow calculation at station 10PT15— if omitted, 10PT15 ■ 1. 

(JLAST Option) 

Stop calculation at station 10PTlb— if omitted, 10rT16 ■ JL. 

(RESTART Option) 

Restart a previously generated case at station I0PT17. 

*N0TEj I0PT9 roust be equal to 3 and KDS must be the same 

value as used in previous run and 10PT14 > 0 In previous run. 
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Card 4 DUCT GEOMETRT CARD F0RMAT (8F10.5) 
10PX3»2 READ DUCT COORDINATES 








Cards 4a DUCT O.D. RADIUS CARDS F0RMAT (8F10.5) 
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Card 4 DUCT GEOMETRY CARO F0RMAT (8F10.5) 

I0PT3»3 STRAIGHT WALL ANNULAR DIFFUSER 
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Card 4 DUCT GE(»1ETRY CARD F0RMAT (8F10.5) 

I0PT3 = 4 COORDINATES ST0RED ON DATA FILE 
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Cards 4a DUCT O.D. RADIUS CARDS F0RMAT (8F10.5) 
10PT3 = 5 CARDS 4a FOLLOW CARD 4 
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Cards 4b DUCT O.D. MIAL STATIONS F0RMAT (8F10.5) 
I0PT3 “ 5 CAkOS 4b FOLLOW CARDS 4a 








Cards 4c DUCT I.D. RADIUS CARDS FjlPMAT <8F10.5) 
I0PT3*2 CARDS 4c FOLLOW CARD 4b 
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INLET/EXIT FLOW DISTRIBUTION F0RMAT (5F10.5) 
I0PT1 » 4 2*KLL CARDS FOLLOWING CARD 5. 
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IF I0PT2 > 0 
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ADD CODE INPUT 
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REFERENCE CONDITIONS F0RMAT (2F10.5,5F6.0,3F10.5) 
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CARD 8 SLOT DATA CARD FORMAT (2110) 
IF (IOPT13.NE.O) 
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CARD 8B DESCRIPTION OF HUB WALL SLOTS FORMAT (510. S) 
NSH CARDS FOLLOW CARD 8A 
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CARD 10 INTERPOLATED OUTPUT DATA CARD FORMAT (IlO, 4F10.5) 





4.3 Output Description for ADO Code 


The printed output on each page of the ADD code Is largely self-explanatory. 

A detailed description of the printed output by page is given together with a sanple 
output . 


Title Page (1) 

Printed by 0UTPUT 

Calculated by duTPUT 

Options None 

Description 

This page presents a list of modifications made to the ADD code together with 
dates and report nximbers. 


1-44 


Echo Print Pages 


Printed by EC0INP 

Calculated by EC0INP 

Options None 

D escription 

The input data is read by subroutine REDINP. This input data is innediately 
printed with input labels by input data card number according to the input data 
sheets given in Section 4.2. This echo print is self explanatory and is inten- 
ded to assist the user in setting up the data cards. 
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Coordinates of Blade Centerline Page (1) 

Printed by FLINE 

Calculated by BLDC0M, FLINE, SLETE 

Options Printed when I0PT2>O 

Description 

The location of the blade centerline in (r,z) coordinates is calculated by 
subroutine BLDG0M. With the (r,z) coordinates known for each blade input data 
point, subroutine FLINE calculates the (n,s) coordinates. Subroutine SLETE 
locates the upstream and downstream blade force calculation surfaces. 


Page 1^ _^£ord^nates_of^ Blad e Center Hne_ 


Heading 

Variable 


Description 



RADIAL LOC. 
OD WALL 

**TCL 

♦ 

Radius-intersection with OD wall 

(ft) 


AXIAL LOC. 
OD WALL 

*TCL 

♦ 

Axial-intersection with OD wall 

(ft) 


WALL DIST. 
on WALL 

’'tcl 

♦ 

Wall distance intersection with 

OD wall 

(ft) 

RADIAL LOC. 
ID WALL 

’^HCL 

♦ 

Radius intersection with ID wall 

(ft) 


AXIAL LOC. 
ID WALL 

*HCL 

♦ 

Axial intersection with ID wall 

(ft) 


WALL DIST. 
ID WALL 

’‘hcl 

t 

Wall distance intersection with 

ID wall 

(ft) 
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Input Flow Data Check Pages (3) 


Printed by WRTCKI 

Calculated by CKINPT 

Options printed when: I0TP1 • 4 or 9 

I0PT5 - 2 


Description 

Subroutine CKINPT checks the Input data used to set up the Inlet and exit flow 
field when I0PT1 ■ 4 or 9 and checks the data used to calculate the blade force when 
10PT5 ■ 2. This subroutine solves the normal momentum equation using the Input data 
to establish radial equilibrium. If the weight flow Is not specified on the Input 
data card, the boundary condition Is set by the static pressure on the ID wall when 
I0PT11 ■ 0 and by the static pressure on the OD wall when I0PT11 *1. If the weight 
flow Is specified, the static pressure is set by the weight flow. In either case, 
the static pressure shown on these pages is that calculated from the normal momentum 
equations. 


Page 1 


Heading 

Variable 

Description 

Y/YT 

Y/Y_ 

Fractional distances across duct 

TOTAL PRES 


Total Pressure (psfa) 

STATIC PRES 

P 

Static pressure (psfa) 

SWIRL ANG. 

a 

Swirl Angle (deg) 

TOTAl. TEMP 

T 

T 

Total temperature (deg R) 

_Pag£ ^ _Che£k_Inpu t_f£r JW^igh^ 

_Rad^al Equilibrium il^PTl_-_91 _ 

Heading 

Variable 

Description 

Y/YT 

Y/Vm 

Fractional distance across duct 

STRM. VEL. 

0. 

Streamwise velocity (ft/sec) 

STAT.. PRES. 

p 

Static pressure (psfa) 

SWIRL VEL. 


Swirl velocity (ft/sec) 

TOTAL TEMP. 

Tt 

Total temperature (deg R) 
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Page 2 Parameters Computed from Input Data 

Heading 

Variable 

Description 

Y/Yy 

Y/Yt 

Fractional distance 

MACH 

M 

Mach number 

STATIC TEMP. 

T 

Static temperature (deg R) 

TOTAL VEL. 

u 

Velocity, ft /sec 

STEM. VEL. 

“s 

Streamwlse velocity (ft/sec) 

TANG. VEL. 


Tangential velocity (ft/sec) 

ROTOR VEL. 


Rotor velocity (ft/sec) 

RELATIVE VEL. 

■ "'b 

Relative velocity (ft/sec) 

RELATIVE ANG. 

B 

Relative angle (deg) 

WT FL^W FRACTION 

W(Y)/W 

Fractional Weight flow 

Page 3 Work Based 

on_Inpu t_Dat£ 

_- 2 i _ _ _ > 

Heading 

Variable 

Description 

Y/YT 

Y/Yj 

Fractional distance 

WORK IN DATA 

Tt2_Tti 

Total temperature rise (deg R) 

WORK IN ROTOR 
ADIAB. EFF. 

'^b(«^2-“4.1> 

(P /P ) 

' 12 ir 

Rotor work input (deg R) 
Adiabatic efficiency 

Y-1 

^I2^’^I1 -1 

ADIAB. LOSS 


Total pressure loss 
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Input Data Pages (A) 


ORIGINAL PAGIK IS 
OF POOR QUALITY 


Printed by 
Calculated by 
Options 

Description 


WRTINP 

FN0RM, FL0W1N 

Pages are printed according to Input options 


The Input data Is printed and labeled including; selected Input options i mesh 
parameters, reference conditions set by BL^CK DATA and subroutine FN0RM, and 
average Inlet flow conditions set by subroutine FL0WIN. 


Page jL _Run_Tj^t le_ 


Heading 

Variables 

Description 

Options used 

1<<PT(^ 

Input options where * 1, 20 

MESH PARAMETERS 

DDS 


Transverse mesh distortion parameter 

KL 


Number of streamlines 

JL 


Number streamwlse stations coarse gri 

KDS 


Number steps/station fine grid 

INLET FL(/W PARAMETERS 


MSI 

Ml 

Inlet Mach number If specified 

ALPl 


Inlet swirl angle If specified 

DSH 

Displacement thickness ID wall (ft) 

DSJ 

Displacement thickness OD wall (ft) 

ANH 

«H 

Power law ID wall 

ANT 


Power law OD wall 

WFLl 

w 

Weight flow If specified (Ib/sec) 

PERFORMANCE POTNT 

WFL(;i 

V p 

Calculated weight flow (Ib/sec) 

REY 

r^PrUy/ur , 

Reference Reynolds number 

DYNPl 

qi 

Mass average dynamic pressure (psfa) 

MACHl 

M , 

Mass average Mach number 

PRESl 

P , 

Mass average static pressure (psfa) 

ATEMPl 

T , 

Mass average temperature (deg R) 

OMEGZ 


Rotor speed (rpm) 

MACHA 

- A 

Area average Mach number 

REYH 

phu/u , 

Reynolds number based on mass average 
flow and inlet height 

B1 

I'l 

Inlet blockage 
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<^fflQlNAL 

POOR 


PAOe IS 

quality 


REFERENCE CONDITIONS 


PRESR 

‘’r 

• 

Reference pressure (psfa) 

TEMPR 


• 

Reference temperature (deg R) 

RHOR 


t 

Kcfi^r*?nco denalty (nlug/it^) 

CP 

‘•P 

t 

Specii Ic heat (ft^/deg) 

CV 

Cv 

1 

Specific heat (ft^/deg) 

VI SCR 

»'r 

• 

Reference viscoHlty (nJug/ft^) 

USK 

Ur 

» 

Reference velocity (ft/»ee) 

KADR 


• 

Reference radiuH (ft) 

SNUK 


• 

Reference hoecd of ttound (ft /sec) 

PR 

I’r 

t 

Prandtl number 

TURBULENCE PARAMETERS 




AKAPPA 

K 

» 

Von Karman conatani 

ACHl 

X 

> 

Clauser constant 

APLUS 

A+ 

1 

Van Driest constant 

PRT 

rj 

t 

Turbulent Prandtl number 

Page 2 Inlet Flow Data If 

ToptI » 

4,8 

Heading 

Variable 

Description 

SPANWISE LOCATION 

Y/Y’ 

T 

t 

Fractional distance 

TOTAL PRESSURE 

**T 

f 

Total pressure (psfa) 

STATIC PRESSURE 

P 

• 

Static Pressure (psfa) 

SWIRL ANGLE 

a 

f 

Swirl angle (deg) 

TOTAL TEMPERATURE 

T 

9 

Total temperature (deg R) 


j0PTl~» 

9 

Heading 

Variable 

Description 

SPANWISE LOCATION 

Y/Y^ 

» 

Fractional distance 

STREAMWISE VELOCITY 

“s 

9 

Streamwise velocity (ft/sec) 

STATIC PRESSURE 

P 

9 

Static pressure (psfa) 

SWIRL VEL. 


9 

Tangential velocity (ft/sec) 

TOTAL TEMPERATURE 

% 

9 

Total temperature (deg R) 
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I’ajjc 3 Strut Da'a _ __ 

Heading Variable Uescriptlun 


ROLOl 

*‘CU) 

t 

Radial location of strut (ft) 

ZCLOl 

*CLO 

• 

Axial location ol strut (It) 

THCLl 

^CLO 

f 

Rotation of strut (deg) 

0M1X;Z1 

iJ 

• 

Rotor speed (rptn) 

NB 

N 

» 

Number of blades 

YCL 

y 

t 

Stacking line Y coordinate of (1 

ALPS 


> 

Stagger angle (deg) 

CH0RD 

c 


Cljord (ft) 

THICK/CHORP 

t/c 

» 

Thickness/chord ratio 

CAMBER 

<^e 

« 

Camber angle (deg) 

xi:l 

X 

« 

X coordinate of tj_ (It) 

Y 

Y/Y-r 

f 

Fractional distance 

BETA!* 

B* 

1 

♦ 

Inlet metal angle (deg) 

BETA2* 

B* 

•> 

• 

Exit metal angle (deg) 


I’age A Strut FL’W Variables 10TP5 ■= 2 

This page is the same as the Inlet Flew Page 2. 
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Duct Geometry Pages (Ti 


Printed by 
Calculated by 
Options 

Description 


WRTGDC 

GDUCT» C00RST, SM00TH 
None 


The output from GDUCT Is printed on page 1 which contains the calculated duct 
coordinates for the ID and OD walls. When I0PT3 ■ 2 or 5 and JLPTS JL» this page 
contains the output from the least squares cubic spline smoothing routine subroutine 
SM00TH. The output from C00RST Is printed on page 2 which contains a shortened 
summary of the calculated coordinates. Note that the wall coordinates for a given 
station number do not agree from page 1 to page 2. On page 1, the coordinates are 
given for equal axial stations when I0PT3 5 and for equal wall arc length when 
I0PT3 “5. On page 2, the coordinates are given for equal stations AS In the com- 
putational plane. The complete set of coordinate data is stored on Unit 9. 


^a_^e_l_ ^alcul^t e_d_I)jw jSemn^tj^ 

Heading Variable 

DUCTI(N) 

STRM. STA 
AXIAL DIST. 

RADIAL DIST. 



Description 

Input parameters 
Streamwlse station ntimber 
Axial corrdlnate (ft) 
Radial coordinate (ft) 


Cal^culat£d_Duc^ Coordinates^ 2 


He«.Jlng 


Variable 


Description 


STRM. STA. 

J 

AXIAL DIST. 

2 

RADIAL DIST. 

r 

WALL DIST. 

X 

CURV. 

k 

1/MET COEF. 

1/h 


Streamwlse station ntimber 
Axial distance (ft) 

Radial distance (ft) 

Wall arc length (ft) 
Curvature (1/ft) 

1/Metrlc coefficient (ft) 


Gap Average Invlscld Flow Page 


Printed by WRTCAL 

Calculated by CALINV 

Options Printed for every JJ station JJ ■ I0PT15, I0PT16 depending 

on print option I0TP4 

Description 

Subroutine CALINV calculates the solution for the approximate inviscid 
rotational swirling flow field and the solution is stored on Unit 22, 


Pa£e_l_ Ga£ Average lnvi£cld_Flow 

Heading Variable Description 


JJ 

JJ 

ZH 

*H 

ZT 

Zj 

Y/YT 

Y/Y. 

TOTAL PRES. 

Pt 

STATIC PRESS. 

P 

SWIRL ANGLE 

a 

TOTAL TEMP. 

Tij 

MACH 

M 

STATIC TEMP, 

T 

TOTAL VEL. 

u 

STEM. VEL. 


TANG. VEL. 


ROTOR VEL. 

T 

Vb 

RELATIVE VEL. 

UA-Vg 

RELATIVE ANGLE 

T 

S 

NORM. VEL. 

«n 


Coarse grid station ntnnber 
Axial location ID wall (ft) 
Axial location OD wall (ft) 
Fractional distance 
Total pressure (psfa) 

Static pressure (psfa) 

Swirl angle (deg) 

Total temperature (deg. R) 
Mach number 

Static temperature (deg, R) 
Total velocity (ft/sec) 
Streamwlse velocity (ft/sec) 
Tangential v.elocity (ft/sec) 
Rotor velocity (ft/sec) 
Relative Velocity (ft/sec) 
Relative angle (deg) 

Normal velocity (deg) 


I-5A 


Wall Bleed Conditions Page (1) 


Printed by WBLEED 

Calculated by WBLEQ) 

Options printed when I0PT8 > 0 

Description 

The wall bleed rate Is estimated a'^prlor:' from the plenum conditions and 
the Invlscid static pressure distribution using subroutine WBLEED. 

Page 1 Wall Bleed Conditions 


Heading 

Variable 

DISCHARGE 

COEFFICIENT 

^IS 

RATIO HOLE AREA 
TO SURFACE AREA 

%S 

PLENUM STAGNATION 
PRESSURE 

Ptp 

PLENUM STAGNATION 
TEMPERATURE 

Ttp 

WALL DISTANCE START 
BLEED 

^L 

WALL DISTANCE 
HUB WALL 


MASS BLEED 
HUB WALL 

^HB 

TOTAL BLEED 
HUB WALL 

'^HB 


Description 

, Discharge coefficient of holes 
f Ratio of hole to surface area 
» Plenum total pressure (psfa) 

, Wall distance start, bleed (ft) 

, Wall distance stop bleed (ft) 
p Wall distance ID (ft) 

, Wall bleed rate ID (Ib/sec/ft^) 
, Integrated bleed (Ib/sec) 
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Heading 

Variable 

HALL DIST 
TIP WALL 


HASS BLEED 
TIP WALL 

«TB 

TOTAL BLEED 
TIP WALL 

^TB 


Description 
Wall distance OD (ft) 

2 

Wall bleed rate OD (Ib/sec/ft ) 
Integrated bleed (Ib/aec) 
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Gap Avera ge Plow Properties Paffaa 

Printed by WRTS0V 

Calculated by S0LVI, TURB, TURB2Q, AMF0R 

Options Printed for every JJ station JJ-I0PT15, I0PT16 depending on 

I0TP4 

I0PT4 > 0 print only page 1 
I0PT4 < 0 print pages 1 through 8 

Description 

Subroutine S0LVI solves the equation of notion for turbulent compressible 
flow using subroutine TURB to calculate the eddy viscosity using algebraic turbulence 
models or subroutine TURB2Q using the two equation turbulence model. Subroutine 
S0LVI also integrates the work Input by the blades and the entropy rise due to 
the dissipation function which is printed on page 6. The solution printed on 
pages 1 and 2 are stored on Unit 8. 




Heading 

Variable 

JJ 

JJ 

JKDS 

JKDS 

AXIAL LOG. 

z 

RADIAL LOG. 

r 

WALL DIST. 

X 

WALL TEMP. 

TV 

WALL BLEED 

»B 

STREAMLINE NO. 

K 

Y/Ytip 

Y/Y» 

STRM. VEL. 

«s 

TANG. VEL. 

NORM. VEL. 

TOTAL VEL. 

u 

SWIRL ANGLE 

a 

MAGH NO. 

M 

TOTAL TEMP. 

Tt 

TOTAL PRES. 


YPLUSH 


YPLUST 

A 


Description 

Streamwise coarse grid number 
Fine grid station number 
Axial coordinate (ft) 

Radial coordinate (ft) 

Wall arc length (ft) 

Wall temperature if specified (deg R) 
Wall bleed (Ib/ft^/sec) 

Transverse grid number 
Fraction distance 
Streamwise velocity (ft/sec) 

Tangential velocity (ft/sec) 

Normal velocity (ft/sec) 

Total velocity (ft/sec) 

Swirl angle (deg) 

Mach number 

Total temperature (deg R) 

Total pressure (psfa) 

Universal distance to first grid point 
ID wall 

Universal distance to first grid point 
OD wall 
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Pa^e_2_ Ga£ Av y^e Fl«#_Pr ©per t^e£ 


Heading 

Variable 

Description 

JJ 

JJ , 

Streamwlse coarse grid no. 

JKDS 

JKDS » 

Streanwlse fine grid no. 

STREAMLINE NO. 

K » 

Transverse grid no. 

Y/YTIP 

Y/Yt 

Fractional distance 

WT. FLOW 

w(Y) 

Weight flow (Ib/sec) 

STATIC PRES. 

P , 

Static pressure (psfa) 

STATIC TEMP. 

T , 

Static temperature (deg R) 

DENSITY 

P » 

Density (slugs/ft^) 

ENTROPY 

HR 

Entropy/gas constant 

STRM. STRESS 


Streamwlse stress (Ib/ft^) 

TANG. STRESS 

4 * 

Tangential stress (Ib/ft^) 

HEAT FLUX 

Heat flxix (Ib/ft/sec) 

Page 3/4 Bounda^ Properties Bub/Tlp Wall 

Heading 

Variable 

Description 

WALL VISCOSITY 


Viscosity at wall (slug/ft/sec) 

WALL HEAT CON. 

♦ 

Conductivity at wall (Ib/sec/deg) 

WALL TEMP 

Tw 

Wall temperature (deg R) 

WALL DENS. 

2/? 2)1/2 ; 

Wall density (slug/ft^) 

WALL STRESS 

Wall stress (Ib/ft^) 

USTAR 

* ^ 

u ^ 

Friction velocity (ft/sec) 

QWALL 

qy(u*)^ » 

Normalizing factor on heat flux 
(Ib/ft/sec) 

STREAMLINE NO. 

K » 

Streamline number 

Y 

y » 

Distance from wall (ft) 

Vise. /WALL Vise 

» 

Effective viscosity ratio 

HEAT CON. 

^E^^w » 

Effective conductivity ratio 

YPLUS 

Y”^ - P„u*y/u„ , 

Universal distance 

UPLUS 

« u/u* » 

Universal velocity 

TPLUS 

1+ - t/T^ 

Universal stress 

QPLUS 

a* - <i/(o„u*’) . 

Universal heat flux 
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Pa&e 5 Gap Average Turbulence Properties 

Heading 

Variable 


Description 

STEAMLIME NO. 

R 

• 

Streamline number 

Y/YTIP 

Y/Yj 

f 

Fractional distance 

KEY. STRESS 


f 

2 2 

Reynolds stress (ft /sec ) 

TURB. K.E. 

u’u’/2 

> 

2 2 

Turbulence kinetic energy (ft /sec ) 

DISSIPATION 

» (tr)^ 


2 3 

Turbulence dissipation (ft /sec ) 

PRANDTL 
MIXING LENGTH 



Prandtl mixing length (ft) 

RICH. NO. 


> 

Richardson number for streamline curva- 

S 

8 


ture 

RICH NO. PHI 


f 

Richardson number for swirling flow 

AvCTage Workj_ Loss^^ 

Efficiency 

Heading 

Variable 


Description 

STREAMLINE NO. 

K 

» 

Streamline number 

Y/YTIP 

Y/Yj 

» 

Fractional distance 

WT. FL0W 
FRACTION 

w(Y) /w 

» 

Weight flow fraction 

ROTOR VEL. 

^B 


Rotor velocity (ft/sec) 

ABS. ANG. 

a • tan“l(u./ug) 

• 

Absolute flow angle (deg) 

REL. ANG. 

a - tan"^(u^-VB>/(u8) 

Relative flow angle (deg) 

WORK IN. 

^T2^^Ti“^ 

• 

Work input 

WORK OUT 



Work output 

LOSS 

1-exp (-AI) 

Y-1 

> 

Loss 

ADIABATIC 

EFFICIENCY 

• 

Adiabatic efficiency 
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truncation 

ERROR 

T 

t 

EROTH 

6(pT) 

Ff 

• 

EPRES 

6P 

r 

f 


6(pus) 


ER0OS 

55s 

9 




EUSU8 

- 2 

^ 

- 2 

9 

EUPUP 

9 




EENTP 

11 

1 

9 


Mss9 average total temperature (deg R) 
Error in PT 

Error in P 

Error in pUg 

, 2 
Error in 

2 

Error in 
Error in I 
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Stream Thrust AvetMe P*ge (1) 


oriqinal 

OF POOR 


page ts 

Q^AUVf 


Printed by 
Calculated by 
Option 


BMF0R 

BMF0R, AMF0R 
None 


Description 

Th. .tre«. thrg.t ««.*. qwntltle. .r. c.lcul.t«l ky Bubroutln. »WR vhlch 
c«vutM « ««.*. M.ch «u»k« «M«h ..tUfU. tht one dlmenelonel energy, 
continuity, and momentum (stream thrust) equations e 


Page 1 Stream Thrust Average Properties 


Heading 

Variable 

Description 

STATION 

J 

Streamwlse station 

NO. 



ZL0C 

Zj^.(Zh+Zi)/2 

Mean axial distance (H) 

AM 

$1 

Average Mach number 

BPR 

>\i 

P 

Average static pressure (psf) 

BPRO 

<V 

Pt 

Average total pressure (psf) 

STRT 

<v 

T 

Stream thrust (lb) 

ARA 

a 

Crossectlonal area (ft^) 

WC 

W 

c 

Choked weight flow (Ib/sec) 

PTL0SS 

(Pti'-Pt)/^T1 

Total pressure loss 

MAMIX 

1 r iT-fidj, 
mJ f 

Total temperature mixing 
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Output Swanarv Pagee (31 

Pr.^kfited by WRTSUM 

Calculated by S0LV1, FAVER 

Option None 

Description 

The mass flow weighted average flow properties are calculated by subroutine 
FAVER and printed on page 1. The remaining flow parameters printed on pages 2 and 3 
are calculated by subroutine S0LVI. 


ORKUNM. PAQC e 
OF POORQUALfTY 


^ 8 ® 


Heading 

Variable 

STEM. STA 

JJ 

AXIAL DIST. 

z * (zg +■ z^)/2 

STRM. VEL. 

“s 

TAN. VEL. 


STATIC PRES. 

p 

STATIC TBIP. 

T 

DENSITY 

p 

MACH No. 

L- 

TOTAL PRES. 

£t 

TOTAL TEMP. 

Tt 

LOSS 

l-exp(Ij^-l) 


w 

BLOCK 

® * ^ ‘ gA(pU), 


Description 

Streamwlse station number 
Axial distance (ft) 
Streamwlse velocity (ft/sec) 
Tangential velocity (ft/sec) 
Static pressure (psfa) 

Static temperature (®R) 
Density (slug/ft ) 

Mach nund>er 

Total pressure (psfa) 

Total temperature (deg. R) 
Loss 

t 

Blockage 


Pa£e_2_ WaUJPjressure md Fr^c^lou_Coeff jx^ent 


Heading 

Variable 


STRM. STA. 

JJ 

» 

AXIAL DIST. 

z 

0 

WALL DIST. 

X 

• 

PRES. COEF. 

(P-Pi)/q. 

> 

STREM. FRICT. 
COEF. 


0 

TANG. FRICT. 
COEF. 


0 

DYNP. PRES. 

Qoo 

0 


Description 

Streamwlse station 
Axial distance (ft) 

Wall arc length (ft) 

Pressure coefficient 

Streamwlse friction coefficient 

Tangential friction coefficient 
Maximum dynamic pressure (psfa) 
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^Convective HeatJTjranaJ^er 


Hefldlng 

Variable 

STRM. STA 

JJ 

AXIAL DIST. 

z 

WALL DIST. 

X 

WALL TEMP. 

T 

LOCAL QW 


TOTAL QT 

q j ^ 


Description 

Streanvlse station 
Axial distance (ft) 

Wall arc length (ft) 

Wall Temperature (deg R) 

Wall heat flux (Ib/ft/aec) 
Integrated heat flux (ft Ib/sec) 


\ 
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1 
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Blade Force Pagaa 


Printed by 
Calculated by 
Options 


F0RC(2) 

F0RC(1), FINVIS, CA8C, (SLADE 
Printed when I0PT2 > 0 


Deecr lotion 

The flow field variables calculated by subroutine FINVIS are printed on 
pages (1) and (2). These variables are printed for both the flov just upstresm 
and just downstream of the blades. The blede row (cascade geometry) used by sub- 
routine CASC Is printed on page 3. The blade force* work and adiabatic efficiency 
calculated by subroutine f0rc are printed on page 4. 


Page 1 Blade Force 

Flow_Vav_labl£s_ 



Heading 
STREAMLINE NO. 

Variable 

K 

0 

Description 
Transverse grid number 


Y/YTIP 

Y/Y^ 

» 

Fractional distance 


STEM. VEL. 

«8 

0 

Streamwlse velocity (ft/sec) 


TANG. VEL. 

u 

0 

Abs' ute tangential velocity 

(ft/sec) 

REL. VEL. 

0 

Relative tangential velocity 

(ft/sec) 

TOTAL VEL. 

0 

Total absolute velocity (ft/sec) 

ABS. ANG. 

a 

0 

Absolute swirl angle (deg) 


MACH NO. 

M 

0 

Absolute Mach number 


TOTAL TEMP. 

Tt 

0 

Total temperature (deg R) 


TOTAL PRES. 

Pt 

0 

Total pressure (psfa) 



Bl®d£ Fprc£ Vajr lab le£ 


Heading 

Variable 


Description 

STREAMLINE NO. 

K 

0 

Transverse grid number < 

Y/YTIP 

Y/Yt 

0 

Fractional distance 

RADIHS 

r 

0 

Radius (ft) 

STAT. PRES. 

P 

0 

Static pressure (psfa) 

STAT. TEMP. 

T 

0 

Static temperature (deg R) 

DENSITY 

P 

0 

Density (slug/ft^) 

ENTR0P/GAS C0NST. 

llfi 

0 

Entropy 

ABS. ANG. 

a 

0 

Absolute flow angle (deg) 

REL. ANG. 

S 

0 

Relative flow angle (deg) 

BLADE VEL. 

''b 

0 

Rotor velocity (ft/sec) 
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Pag 0 3 Cascade 

Paranecaro 



Heading 

Variable 


Description 

STREAMLINE NO. 

K 

0 

Transverse grid number 

RADIUS 

r 

• 

Radius (ft) 

THICK/CH0RD 

t/c 

f 

Thickness/chord 

GAP 

s 

0 

Gap (ft) 

S0LD 

0 

0 

Solidity c/g 

CHORD 

C 

0 

Chord (ft) 

CAMBER ANC. 

♦c 

0 

Circular arc camber angle (deg) 

STAGGER ANG. 


0 

Stagger angle to axis (deg) 

LOSS COEF. 

^B 

0 

Loss coefficient 

Pa&e 4 Blade Force* Work* Efficiency 


Heading 

Variable 


Description 

STREAMLINE NO. 

K 

9 

Transverse grid number 

LIFT C0EF. 

^L 

0 

Lift coefficient 

DRAG C0EF. 

Cd 

0 

Drag coefficient 

STRM. F0RCE 

fs 

0 

StreaiDwise force/span, Ib/ft 

TANG. F0RCE 


9 

Tangential force/span* Ib/ft 

WORK INPUT 

^p^^T2“^Tl^ 

0 

o 2 

Work input* ft*/ sec 

TT2/TT1-1 

Tt2/T^1-1 

0 

Total temperature increase 

PT2/PT1-1 

P„,/P„ -1 
T2 T1 

9 

Total pressure increase 

ADIAB. EFF. 


9 

Adiabatic efficiency 


PAGE IS 
OK HiOR QUALITY 
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rrlnted by BLPARM 

Calculated by BLPARM 

options 

ivscrijvtion compared with the viscous 

The invlscid flow field solution boundary layer is determined using a 

solution stored on Unit 8 and the "^^^"^^^.'’^omenrum thicknesses are calculated 
vortlcitv criteria. Then the displacement an swirling flow with 

using the definition the definitions have an extended 

thick boundary layers and normal pressure gr 

form. 


Page_l/2_ Oi 

Heading 

Z 

DEL 

MAai 

DELS 

THET 

Hi: 

RETH 

DELH 

THEP 

CF 


5rm.n5io»al B 0 undary_.Uyer_Paramy.r 5 . 


Variable 


Description 

Axial distance (ft) 

Boundary layer thickness (ft) 

Edge Mach number 
Displacement thickness (ft) 
Momentum thickness (ft) 

Shape factor 

Reynolds number based on 6 
Swirl displacement thickness ( 1 

Swirl momentum thickness (ft) 
Skin friction coefficient 



Interpolated Output data 

Printed by 
Calculated by 
Options 


pages (7) 

FLINE, WRTF0U 
FLINE, F0UTP 
Printed when NLAST > 0 


Description 

The (n,s) coordinates for the output data line are calculated and v 

...brouttne PUNE. The lnt.rpoUt«i solution on the output d.t. fln'^. 

SillW sh»rth?”,‘« "'.‘"‘1 Print«i outfit JtL *" 

in three di«went coordinate systems. Pages 2 and 3 show the solution in the 
(n.s.j) or add c^e coordinate system. Pages 4 and 5 show the solution in the 
th* *fv coordinate system. Pages 6 and 7 show the solution In the 

the (Y.X.«) coordinate system where Y is tangentent to the output data line and 
X Is normal to the output data line. 


Page 1 _Coord^nates_oX OalPH. t_Da ta Line_ 

Variable Description 


RADIAL L0C. 
OD WALL 


Radius Intersection with OD wall (ft) 


AXIAL L0C. 
OD WALL 


Axial intersection with OD wall (ft) 


WALL DIST. 
OD WALL 


*WT 


» OD wall distance to intersection (ft) 
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■p”gg 2 ""cap Average Flow Properties at Output Data Station 


(N,S)_Coordif«te Systan_ 

Heading Variable 

AXIAL L0C. Zg 

HUB 


RADIAL LOG. 


HUB 

WALL DIST. 

% 

HUB 

♦WALL TEMP. 

% 

HUB 

WALL BLEED 


HUB 

AXIAL L0C. 

h 

TIP 

RADIAL L0C, 


TIP 

WALL DIST. 


TIP 

♦WALL TEMP 

Tt 

TIP 

WALL BLEED 

0 

K 

T 

♦Note: For adiabatic 

vails 

STREAMLINE NO. 

L 

Y/YTIP 

Y/Y, 

STRM. VEL. 


TANG. VEL. 


N0RM. VEL. 

“n 


Description 

, Axial wall coordinate ID (ft) 

, Radial wall coordinate ID (ft) 

, Wall distance ID (ft) 

, Wall temperature ID (deg. R) 

, Wall bleed ID (Ibn/ft^/sec) 

, Axial wall coordinate CD (ft) 

, Radial wall coordinate OD (ft) 
, Wall distance OD (ft) 

, Wall temperature OD (deg. R) 

, Wall bleed (» (Ibn/ft^/sec) 

0 . 

, Output data point number 
, Fractional distance 
, Streamwlse velocity (ft/sec) 

, Tangential velocity (ft/sec) 
Normal velocity (ft/sec) 
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Heading 

Variable 


Description 

RADIAL L0C. 
ID HALL 

’^WH 

• 

Radius intersection with ID wall (ft) 

AXIAL L0C 
ID WALL 

*HH 

» 

Axial intersection with ID wall (ft) 

WALL DIST 
ID WALL 

*HH 

» 

ID wall distance to intersection (ft) 

P0INT NO. 
L 

L 

» 

Output data point ninber 

RADIAL L0C. 
FT. 

'l 

P 

Radial location of point L (ft) 

AXIAL L0C 
FT 

*L 

P 

Axial location of point L (ft) 

STRM. DIST. 
FT 

*L 

P 

Streamwise distance to point L (ft) 

STRM. C00R. 
S 

®L 

P 

Streanwlse coordinate of point L 

NORM. C00R 
N 

“l 

P 

Momal coordinate of point L 

*STRM. STAT. 

j 


P 

Streamwise mask point nusiber 

*N0RM. STAT 

K(sL,nL) 

P 

Normal mask point number 


K 

*The point lies between (J» J+1) and (k^krt-l) • 
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Heading 

Variable 

Description 

TOTAL VEL. 

u 

» Total velocity (ft/sec) 

SWIRL ANGLE 

o»tan"l(UA/Oa) 

, Swirl angle (deg.) 

MACH NO. 

M 

» Mach number 

TOTAL TEMP. 

Tt 

» Total temperature (deg. 

TOTAL PRES. 

Pt 

, Total pressure (psta) 


Page 3 Gap Average Flow Properties on Output Data Line 
(n,£)_Coord^nate ^ystra 


Heading 

Variable 

STRMLINE NO. 

L 

Y/YTIP 

Y/Y» 

WT. FLOW 

W 

STATIC PRES. 

P 

STATIC TEMP. 

T 

DENSITY 

P 

ENTROPY 

I/R 

STRM. STRESS. 

ans 

TANG. STRESS 

®n^ 

HEAT FLUX 

q 


Description 

, Output da^«> point no. 

, Fractional distance 
, Weight flow to point L (Ib/sec) 
, Static pressure (Ib/ft^) 

, Static temperature (deg. R) 

, Density (slugs/ft^) 

, Entropy (dimensionless) 

, Streamwise stress (Ib/ft^) 

, Tangential stress (Ib/ft^) 

, Heat flux (Ib/ft/sec) 
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Page 4 Gap Average fIom Properties at Output Data Station 
(r,£)_Coord^nat£ 


Heading 

Variable 

STREAMLINE NO. 

L 

Y/YTIP 

Y/Y^ 

AXIAL VEL. 


TANG. VEL. 


RADIAL VEL. 

“r 

TOTAL VEL. 

u 

TAN(A)»UP/VZ 

a“tan“^(U,j,/U2) 

MACH NO. 

M 

TOTAL TEMP. 

Tt 

TOTAL PRES. 

Pt 


Page 5 

Gap Average Flow Properties 
(^z) jC£Ojrdinat£ Sy£tCT_ _ _ 

Heading 

Variable 


Description 

, Output data point no. 

» Fractional distance 
, Axial velocity (ft/sec) 

, Tangential velocity (ft/sec) 
, Radial velocity (ft/sec) 

, Total velocity (ft/sec) 

, Swirl angle (deg.) 

, Mach number 

, Total temperature (deg. R) 

, Total pressure (psfa) 


at Output Data Station 


Description 


STREAMLINE MO. 
ANGLE T 
WT FLOW 
STATIC PRES. 
STATIC TEMP 
DENSITY 
ENTR0PY 
SIGRF 
SIGRP 
HEAT FLUX 




q 


, Output data point no. 

, Angle between S and Z coordinate (deg.) 
, Weight flow (Ib/sec) 

, Static pressure (psfa) 

, Static temperature (deg. R) 

, Density (slug/ft^) 

, Entropy (dimensionless) 

, Stress component (psfa) 

, Stress component (psfa) 

, Heat flux (Ib/ft/sec) 
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Page 6 Gap Average Flow Propertlea at (Xitput Data Station 
( 3^^ 3^) ^Coord limt e Syatm^ 


Heading 

Variable 


Description 

STREAMLINE NO. 

L 

» 

Output data point no. 

Y/YTIP 

Y/Yj 

f 

Fractional distance 

UX VEL. 


f 

velocity (ft/aec.) 

TANG. VEL. 

UA 

# 

Tangential velocity (ft/aec.) 

UY VEL. 

4 

♦ 

Uy velocity (ft/sec.) 

TOTAL VEL. 

u , 

f 

Total velocity (ft/sec) 

TAN(A)« UP/DX 

o»tan"^(®a/^t) 

f 

Swirl angle (deg.) 

MACH NO. 

M 

• 

Mach number 

TOTAL TEMP. 

^T 

f 

Total temperature (deg. R) 

TOTAL PRES. 

Ft 

> 

Total pressure (psta) 


Page 7 Gap Average Flow Propertlea at Output Data Station 


(y,x) Coordinate System 
Heading Variable 


Description 

STREAMLINE NO. 

L 

» 

Output data point no. 

ANGLE T 

0 

f 

Angle between X and Z coordinates (deg.) 

WT. FLOW 

W 

f 

Weight flow to point L (Ib/sec) 

STATIC PRES. 

p 

» 

Static pressure (psfa) 

STATIC TEMP. 

T 


Static temperature (deg. R) 

DENSITY 

P 

f 

Density (deg/ft^) 

ENTROPY 

I/R 

$ 

Entropy (dimensionless) 

SIGXY 

V 

> 

G^y stress (psf) 

SIGYP 

°y^ 

f 

Gy^ stress (psf) 

HEAT FLUX 

q 

f 

Heat flux (Ib/ft/sec.) 
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4.4 Diagnoetlcs for ADD Code 

Numerous checks are made during the course of the calculation. If a minor 
error occurs, a DIAGN0STIC message is printed and the calculation continues. 

If a fatal error occurs, a DIAGN0ST1C message is printed and the calculation stops. 
A description of the DIAGN0STICS is given in this section. The DIACN0STIC message 
is always in the form: 

**DIAGNOSTIC NO. XX FOR ANNULAR DIFFUSER DECK** 

where xx refers to one of the errors listed. It should ha noted that numerical 
values printed with the DIAGNOSTIC massage will ha in dimenslonlass form or in 
English units. 

1) I0PT3 OUTSIDE RAIKSE OF ALLOWABLE DUCT OPTIONS 

This error is detected in Subroutine ALTMN. The input option snist be between 

1 < 10PT3 £ 5. 

2) No solution exists in AMFOR 

This error is detected in Subroutine AMFOR. This subroutine solves the Mach 
nvimber function 

for M given N. The function has a maximum at M • 1. Hence 

- 1/2 

N (l)« [2(l + y)] 

Solutions do not exit for values of N > N(l). 


3) MASS FLOW EXCEEDS THE MAXIMUM MASS FLOW POSSIBLE 

This error is detected In Subroutine AMINLT which solves the Mach number 
function y» I 

/ ^ 

.. A .. /. . _K. L ..2 1 


N*M + M*) 


for M given N. This function has a maximum for M > 1 given by 

r*l 

N(l) 

corresponding to choked flow. 


4) I SHAPE AND 10PT2 ARE NOT CONSISTENT 


This error is detected in subroutine CASC. For blade and strut calculations 
ust? only I0PT2-3 with any ISHAPE, where 
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3 < ISHAPE s 6 
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Otherwise, the calculation will stop. 


5) FOR BEST RESULTS ADD A STRAIGHT ANNUUR CHANNEL INLET 


This error is detected In Subroutine Cddni t ..u 

duct coordinates, it Is assumed that the lnleJ”h construction of the 

..ccon .0 ... 



6) prooium assumes inlet flow has curvature 

This error Is detected In Subroutine CWRL. Same 


/ 5 > 





as diagnostic 5. 


7) WALL CURVATURE IS TOO LARGE AT STATION X. 

This error is detected in Subroutine C00R1 usually if the duct has a 
discontinuous change in wall curvature such as shown in the figure below. 



Discontinuous Change in Wall Curvature 


8) Not Used 

?) GREATER THAN 1.0 PERCENT NORMAL PRESSURE GRADIENT ERROR RECALCULATE STATIC 
PRESSURE 


This error is detected in Subroutine ERPIN. This subroutine integrates 
the radial tnomentum equilibrium equation. 





L 

R dn 



and compares (Pj - Py) to that computed for the input inlet flow (P-p - Pjj)j 
If the error given by 




I - 


tP,- P«l 


I 


is greater than 0.01. the input initial static pressure distribution is 
replaced by the above pressure equation and the inlet flow is recalculated. 
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10) Not Used 


11 ) MASS FLOW REQURED EXCEEDS MAXIMUM MASS FLOW POSSIBLE 


This error is detected in Subroutine CKINPT. If it is determinud that choked 
flew exists in the duct, this diagnostic will be printed; the weight flow must 
be reduced. 

12) PRESSURE RISE EXCEEDS PERMISSIBLE PRESSURE RISE 

This error is detected in Subroutine CKINPT snd indicates that the deck cannot 
calculate properly the initial flow profiles. Check input for errors. 

13) ITERATION OF BACK PRESSURE CALCULATION FAILS TO CONVERCE 
This error is detected in Subroutine FINVIS. 

In the calculation of strut forces, it has been assumed that the strut exit 
flow is subsonic and unseparated (l.e., Ug >0). If these conditions are 
violated, no solution can be obtained. The calculation will stop. 

14) BOUNDARY LAYER TOO THIN FOR MESH SPACING 

This error is detected in Subroutine FL0WIN. The viscous flow calculation 
requires a finite initial boundary layer thickness. In addition, it requires 
enough stesh points to describe the inlet boundary layer velocity profile. 

The deck assumes arbitrarily that at least five mesh points are required. 

Thus, if this diagnostic occurs, increase the number of mesh points, KL, 
increase the mesh distortion parameter, DDS, or Increase the assumed inlet 
boundary layer thickness. If DDS is input equal to zero, the program auto- 
matically sets the mesh distortion parameter to the appropriate value for 
turbulent flow. 

15) TOTAL PRESSURE IS LESS THAN STATIC PRESSURE 

This error is detected in Subroutine FL^WIN. A check is made on the input 
data for I0PT1> 4 to make sure that P^ > P. 

16) INPUT DATA NOT IN RADIAL EQUILIBRIUM CORRECTIONS APPLIED TO STATIC PRESSURE 

This. error is detected in Subroutine FL0WIN. A check is made of the input 
static pressure data for I0PT1*4. If the static pressure data are not in 
radial equilibrium, it is assumed that the static pressure data are in error 
and that the other inlet data are correct. Then the static pressure profile 
is computed from 
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d n 

d 7f 


2 


JL{zL 

r-i Lxv 


dv 2 
— COS a 


-L Ifi 

XR dn 


sin 



with the ID wall static preaaure as a boundary condition. 


W POOH ijUAltn 


17) INPUT DDS MUST BE SPECIFIED 


This error Is detected In Subroutine FNdRM. At this time there Is no 
algorithm to select automatically the mesh distortion parameter DDS for 
laminar flow. 


18) BLADE DATA ERROR IN CKINPT ROUTINE 

This error Is detected In Subroutine CKINPT. Blade data have been Input 
Incorrectly and must be rearranged with Y Increasing. 

19) NO UNIQUE SOLUTION FROM MINVRT 

This error Is detected In Subroutine MINVRT. If the matrix used to solve 
for the turbulent flow solution is singular, no solution can be obtained. 
This situation may occur due to numerical truncation errors. 

20) LEADING OR TRAILING EDGE INDEX OF STRUT OUT OF RANGE 

This error Is detected in Subroutine SLETE. In order to compute blade 
forces, the strut roust be located entirely within the duct length. This 
problem may be eliminated by extending the duct as shown in the figure. 
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21) SLOT INPUT NOT IN INCREASING ORDER 
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This error is detected In Subroutine SL0TA. 

The slot Input dnta must b«- arranged In order of increasing axial distance. 
Check input data. The calculation stops If this error is detected, 

22) CHOKED FLOW IN SLOT. 

This error is detected in Subroutine SLTFL0, The slot weight flow is 
determined by the ratio of the stagnation pressure of the slot coolant 
fluid to the local wall static pressure. If this pressure ratio is too 
large the flow may be choked at the slot inlet. The calculation will stop. 

23) BOUNDARY LAYER OVERLAP OR TOO LARGE 

This error is" detecte'id in Subroutine FL0W1N, For internal flow, the sum 
of the boundary layer thicknesses on the hub and tip walls must be less 
than the duct inlet height. Check input data. 

24) SET TOTAL TEMPERATURE. PRESSURE, ANGLE TO VALUE AT EDGE OF BOUNDARY LAYER 
- CORRECTION APPLIED 

This error ic detected in subroutine FL0WIN, For I0PT1-4, the calculated 
boundary layer profiles are matched to input inlet flow profiles. 

A good match requires that the stagnation pressure, Pij^, bo constant in the 
experimentally determined boundary layer region. as shown in the figure 
(dashed line) . 



Constructing the Inlet Flow 



25) TRUNCATION ERROR CANNOT BE REDUCED BY STEP SIZE 


This orror 1 b detected In Subroutine S0LVI. If the step size parameter (KDS) 
ia not specified, it is selected automatically by checking the truncation 
error at each step. When an Instability occurs, the program attempts to 
reduce the truncation error by reducing the streamwlse step size. If the 
truncation error cannot be reduced below a mlnumum value, the calculation 
stops. 

26) NUMERICAL INSTABILITY 

This error Is detected In Subroutine FC^^RCT and Subroutine S0LV1 and is an 
indication that the program has calculated negative temperature or pressures. 
The calculation stops if this problem occurs. 

27) RHOeX ITERATION DID NOT CONVERGE, ERR ■ 

This error is detected in Subroutine F0RCE. In determining the blade force, 
an iteration scheme is used to determine the downstream static pressure. If 
this iteration fails to converge, this diagnosti'; is printed together with the 
maximum error found in the iteration. The calculation, however, is not 
terminated . 

28) 10PT3 <= 2 OPTION NOT IN USE 

This error is detected in Subroutine F0RCE but this option has been deleted 
from the current version of the ADD code. 

29) SOLUTION REQUIRES REVERSE FLOW, INCREASE WFLOW 

This error is detected in Subroutine CKINPT. For flows with radial pressure 
gradients, there is a minimum weight flow below which reverse flow exists. 

This problem can be corrected by increasing the weight flow. The calculation 
will stop. 

30) LOAD DOWNSTREAM FLOW DATA CARDS 

This error is detected in subroutine CALINV and Indicates that the downstream 
flow data cards, required by 10PT1»A or 9, have not been entered. The 
calculation will stop. 

31) SOLUTION FOR BLADE FORCE DOES NOT EXIST 

This error is detected in subroutine F0RCE. The blade force cannot be 
calculated because no iuviscid flow solution can be calculated. (Same 
as OIAGN0STIC 29). The calculation will stop. 
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32) GRADIENT OF METRIC COEFFICIENT • 

FOR BETTER RESULTS ADD STRAIGHT CHANNEL INLET 
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lltl in Subroutine C00R4. It is assumed that the Inlet 

JrLa problems, add a straight annular section 

to the inlet. The calculation will continue. 

33) INPUT TOO LARGE FOR COLE’S LAW 
SET N < 

This error is detected In subroutine FL0WIN. 

Cole’s friction law requires a certain relationship Hi, • H 
there 1. upper bouiui of n < 10. For . .olotlon to «ie. “ ‘ 




Setting 


Then 


LS^l 

(cue 

u* 


n <^—-1 

”tnln 

The calculation will stop. 

34) WEIGHT FLOW ITERATION MAY NOT CONVERGE IN SUBROUTINE CKINPT CHECK INPUT DATA. 

error is detected in Subroutine CKINPT. The weight flow iteration that 

flow’is^hl pressure may not converge if the free stream invlscid 

oieLuro ^ ^ distorted. An input flow which is more uniform in stagnation 

pressure is required. The calculation will stop. 


35) WFLI and I0PT11 OPTIONS INCOMPATIBLE. 


to «l8ht flow cannot be specified 

for external flow. The calculation will stop. 


36) I0PT1 = 1 or I0PT1 = 2 OPTIONS NOT USED. 


This error is detected in Subroutine ALIMN. The options I0PT1 - 1 and I0PT1 = 2 
have been deleted from the code. x ana i«»pti = 2 
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37) CONFLICT OF OPTIONS, I^PT14 < 0 IMPLIES SEPARATION AND GLOBAL ITERATIONS. 
AUTOMATIC STEP SIZE ALGORITHM CANNOT BE USED. 

This error Is detected In Subroutine S(2lLVI. VAten performing a global 
Iteration (I^PTIA < 0), the same number of streamwlse steps must be used 
for each Iteration. Hence the automatic step size algorithm for the 
streamvlse direction must not be used. The ralculatlon will stop. 

38) LET TOO SMALL FOR TURBULENT BOUNDARY LAYER 
RET - XXXX 

This error Is detected In Subroutine CFC0LE. Cole's skin friction law Is 
not valid for Reynolds numbers based on momentum thickness. RET less than 
1000. For RET less than 1000, the boundary layer Is laminar. The 
calculation will stop* 

39) DATA LINE DOES NOT INTERSECT WALL WITHIN MESH 

This error Is detected in subroutine ALINE. The first and last points on 
the Input data line must lie outside the duct and second and next to last 
must lie Inside the duct. The calculaticn will stop. 

40) INP'JT DATA LINE DOES NOT INTERSECT WITHIN MESH 

This error is detected in Subroutine CR0SS2. Check points on Input data 
line to make sure that they lie Inside the computational mesh. The 
calculation will stop. 

41) CAMBER ANGLE <{> IS OUTSIDE LIMITS OF CORRELATION 

This error is detected In Subroutine CASC. If ISHAPE ** 5 or 6, then for 
correlations in NACA SP36, > 0. The calculation will continue but the 

user Is outside limits of correlations. 

42) IMSL LIBRARY FAILURE NO. XXX 

This error Is detected in Subroutine »I00TH. Check input coordinates 
specified on card(s) 4. Check ISCVKU routine from IMSL library. 
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4.5 Debug Options for ADD Code 


When set equal to unity • these options allow interaedlate results calculated 
by thp subioutine Indicated to be printed as an aid in debugging s troublesone 
cn<?e. Note that these outputs are not converted to units and reference 
otifir b<> ttade to the source code for interpretation of printout. 


OPTION 

SUBROUTINE 

OBJECTIVE OF SUBROUTINE 

IDB6I 

TURB 

Calculates two-layer turbulence oodel 

I0BG2 

PC0RCT 

Calculates shear stresses and heat fluxes at 
each station 

IDBG3 

FL0WIN 

Generates initial flow profiles 

1DBG4 

SLTFL0 

Calculates slot inlet flows 

IDBG5 

S0LVI 

Calculates viscous flow solutions 

IDBG6 

C00R 

Generates required geonetric parameters 

IDBG7 

F0RCE 

Calculates forces generated by struts and 
blades 

IDBG8 

MINVRT 

Inverts a matrix 

IDBG9 

SM00TH 

Smooths duct contour read via I0PT3 ■ 2 

IDBGIO 

GDUCT 

Calculates duct geometry 

IDBGll 

SLTFL0 

Obtains additional information from SLTFL0- 
see IDBG4 

IDBG12 

S0LVI 

Obtains additional information from S0LVI - 
see IDBG5 

1DBG13 

CKINPT 

Checks inlet flow input for errors 

IDBG14 

S0LVI* 

Debugs the algorithm that automatically 
computes the maximum step size in the stream- 
wise direction while assuring computational 
stability. 

1DBG15 

IDBG16 

Not used. 

Specifies number of streamlines to use 
in COORST calculations (Default 25). 

1DPG17 

Not used. 
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A* 6 Sample Input for ADD Code 


Two aarople inputs to the ADD code are presented on the following pages. 
These cases correspond to the design studies described in Sections A.l and 4.2. 
The first sample is the input for the Axis,Tnmetric Compressible Curvature Case 
and the second sample is the input for the Separated Flow Case. 


Axlsywmetrlc Compressible Curvature Case 


The option card (line 2) Indicates that the inlet flow is to be 
calculated assuming a constant stagnation pressure and stagnation temperature 
in the core flow (I0PT1-3). The duct geometry is to be read from input data 
cards (I0PT3-2). On the mesh parameter card (line 3), the default mesh 
distortion parameter (DDS-0) has been selected but the streamwlse step size 
parameter has bs?n input KDS-5. From line 3. it is noted that the duct 
coordinates at JLPTS»100 equally spaced axial stations are to be read. - e 
length ot the duct is 4.19375 ft (line 4). Llnet- 5 through 17 contain 1 . 
data points for the tip radii and lines 18 through 30 contain 100 data polu.s 
for the hub radii. The inlet Mach number is 0.7 (line 31). 


Separated Flow Case 


The option card (line 2) Indicates that the inlet flow is to be calculated 
assuming a constant stagnation pressure and stagnation temperature in ^ 
flow (10PT1-3). The duct geometry is to be read from input data cardq (W^3 2). 

On the mesh parameter card (line 3), the mesh distortion parameter 
and the streamwlse step size parameter has been input KDS«*2. Fr^ line 3, 
is noted that the duct coordinates at JLPTS-80 equally spaced axial «tatJ^ons 
are to be read and that the least squares spline smoothing routine will be 
used (JLPTS ^ JL), The length of the duct is 1.764 (Itoe 4). Lines ^ throug 
14 contain 80 data points for the tip radii and lines 5 through 24 contain 
data points for the hub radii. The inlet Mach number is 0.287 (line 19), 
stagnLion pressure and the stagnation temperature ^ at^spheric conditions 

(Ft - 1 atm and Tt - 519°R. For this card I0PT14-1 which indicates that 
global iterations will be used and for the first pass, the convection teras will 
L set to zero in regions of separated flow. If I0PT14>1 then windward differencing 
will be used in calculating the convection terms. 
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